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DEGENERATE DIFFUSION ARISING IN TECHNICAL TEXTILE 

INDUSTRY 
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Abstract. We study a stochastic equation modeling the lay-down of fibers in the pro- 
duction process of nonwovens. The equation can be formulated as some manifold-valued 
Stratonovich stochastic differential equation living on R'' x S"*"^, d > 2. Especially, 
^N) we study the long time behaviour of the stochastic process. Demanding mathematical 

T-H difficulties arising due to the degeneracity of the lay-down equation and its associated 

^^ generator. We prove strong mixing properties by making use of the hypoellipticity of the 

Cm generator and a new version of Doob's theorem derived recently in |GN12) . Moreover, 

^ we show convergence to equilibrium exponentially fast with explicitly computable rate of 

convergence. This analytic approach uses powerful modern Hilbert space methods from 
the theory of hypocoercivity developed in [DMSIOJ . Summarizing, we give interesting 
mathematical applications of geometric stochastic analysis to real world problems. 
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1. Introduction 



This article is about the mathematical analysis of new fiber lay-down equations devel- 
oped recently in j KMWlS] . Herein fiber lay-down processes arise in the production process 
of nonwovens and the expression is used for the description of the forms generated by the 
stochastic lay-down of flexible fibers. The understanding, optimization and mathematical 
simulation of such fiber webs is of great industrial interest, we refer to |KMW09] and 
^7^ references therein. Areas of application where these stochastic lay-down processes can be 

P^ observed include e.g. composite materials (filters), textiles, as well as the hygiene indus- 

(^ try. In |MW06] a general mathematical model describing the full fiber spinning process is 

*^ introduced and is nowadays implemented in the software tool FYDIST developed at the 

Fraunhofer ITWM, Kaiserslautern. Nevertheless, the numerical simulation leads to exces- 



^n sively large computation times. Hence it would be desirable to have simplified stochastic 

models at hand simulating a virtual fiber web in a fast and efficient way. The first of such 
surrogate models, the basic two-dimensional model, is developed in |GKMW07j and reads 
^^ as some Ito stochastic differential equation (SDE) in M^ of the form 

"^ dCt = r{at)dt (1.1) 

Here T{a) = (cos(a),sin(a))^, r-*- = 4^, o" is a nonnegative constant and VF is a standard 
one-dimensional Brownian motion. $ : M^ — )■ M is a suitable function called the potential. 
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This basic model has been extended in |KMW12] to the more reaHstic three-dimensional 
case which serves as starting point for all our mathematical studies in the underlying 
article. And exactly this extension to the three-dimensional case requires the usage of a 
differentialgeometric language. 

At first we discuss the geometry underlying the previously mentioned equation intro- 



duced in |KMW12] . Starting with the two-dimensional model (1.1) we present a new 
differential geometric derivation and formulation of the model from |KMW12] in each di- 



mension d G N, d > 2. Therefore we formulate Equation (1.1 ) first in its most natural way 
on M? X M/27rZ. The latter space is diffeomorphic to M^ xW where S stands as usual for 
the unit sphere. Consequently, we get the analogue equation of (1.1) living on M^ x S^. 



The resulting equation can then directly be translated to higher dimensions and yields 
the final d-dimensional fiber lay-down model given as some manifold-valued Stratonovich 
SDE with state space M°' x S^'-^ by 

dCt = vtdt (1.2) 

dvt = -{I -vt(E> vt)V<^{S,t) dt + a{I -vt(E) vt) o dWt- 

Here x (^ y = xy , ^ G C^{M. ) and VF is a standard d-dimensional Brownian motion, 
see Section [3] for details. This precise fomulation is also necessary for all the forthcoming 
analysis and we strongly believe that the geometric derivation can be helpful for every 
applied mathematician who aims to derive similiar manifold-valued stochastic equations. 



Furthermore, we provide concret numerical simulation formulas of (1.2) in specific local 



coordinate systems, so called local SDEs. As consequence, for d = 2 we obtain back (1.1 ) 



and for d = 3 the local SDE reduces to the three-dimensional model derived in |GKMW07 



We finish the geometric discussion in Section |4j Therein we further discuss basic existence 
statements of all occuring stochastic equations by constructing global solutions to the 
basic two-dimensional model (1.1) with state space M'^ and the general geometric fiber 



lay-down model (1.2) living on R x S' 



jd-i 



Of course, the cases d = 2 and d = 3 are the physical relevant ones. Nevertheless, for all 



the forthcoming analysis, it is elegant to study Equation (1.2) in its general mathematical 
form. We emphasize that our upcoming stochastic and functionalanalytic considerations 
of the d-dimensional fiber lay-down model are done in a coordinate free way. 



The density of the stochastic process solving (1.2) satisfes the associated Fokker-Planck 
evolution equation. An important criterion for the quality of the fiber web and the resulting 
nonwoven material is how fast the process converges towards its stationary state. Of 
essential interest is therefore the speed of convergence. Moreover, from a practical point 
of view, process parameters should be adjusted in order to obtain optimal convergence to 
equilbrium. 

This motivates that we devote the main interest in the underlying paper to the study 
of the long-time behaviour of the d-dimensional fiber lay-down equation. Demanding 



mathematical difficulties are then occuring due to the degeneracity of (1.2). The con- 
vergence to equilibrium of the „flat" two-dimensional fiber lay-down model is already 
analyzed in several articles. Basing upon the theory of Dirichlet forms and operator semi- 
groups, the approach of [GKOS^ gives an ergodic theorem and establishes explicit rates 
of convergence. The underlying object of study therein is the two-dimensional fiber lay- 
down Kolmogorov PDE and its corresponding Kolmogorov operator. Another approach 
is presented in |DKMS12] . There the authors analyze the two-dimensional fiber lay-down 
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Fokker-Planck PDE together with its associated generator and prove convergence to equi- 
Hbrium with an exponential, exphcitly computable rate of convergence. This approach 
uses modern methods from |DMS10] in which a new hypocoercivity theory in a Hilbert 
space setting is developed. For a general study of the theory of hypocoercivity, the reader 
may consult | Vil09j . Finally, the quite recently published article |KSW11] uses a prob- 
abilistic approach. The authors are able to derive strong mixing properties under weak 
assumptions on the potential and, assuming some stronger conditions, they can prove 
geometric ergodicity of the two-dimensional fiber lay-down process. 

In Section [5] we start studying the long time behaviour of our general, geometric d- 



dimensional fiber lay-down process by proving strong mixing properties, see Theorem 5.3 
The approach is disjoint and different to the one given in the two-dimensional setting 
in |KSWll] . and moreover, it has even been developed at the same time independently. 
Our approach works for smooth potentials $. This is a slightly stronger assumption 
as made in [ KSWll] in the two-dimensional case. Nevertheless, our strategy applies to 
the full manifold-valued d-dimensional fiber lay-down SDE containing especially the two- 
dimensional situation. We make use of a new version of Doob's theorem, see |GN12) . 
which perfectly fits into the fiber lay-down setting. 

Afterwards, we switch to functional analysis and apply the fascinating and powerful 
hypocoercivity theory from |DMS10) once more, see Section pi We generalize the strategy 
in |DKMS12] to the d-dimensional setting requiring some differential geometric tools. The 
object of interest in this section is the hypoelliptic Kolmogorov operator associated to 



(1.2) given by 

2 

L = vV^-{I-v(g) v)V<^{C) -Vv + Y ^^"-^ 

which is analyzed in an appropriately chosen L^-space. In particular, we make use of 
modern entropy methods. We obtain convergence to equilibrium exponentially fast with 
explicitly computable rate of convergence, see Theorem |6.4[ 



Moreover, in Section |2] we introduce some basic geometric language and give a short 
introduction to the concept of manifold-valued Stratonovich SDEs. Further useful state- 
ments needed for our analysis are proven in the Appendix. 

Finally, the progress achieved in this paper may be summarized by the following list of 
main results: 

• Differential geometric derivation and formulation of the fiber lay-down model in 
each dimension as some manifold-valued Stratonovich SDE. 

• Construction of strong solutions to all occuring fiber lay-down equations. 

• Strong mixing properties of the geometric fiber lay-down process under weak as- 



sumptions on the potential, see Theorem 5.3 by making use of a new version of 
Doob's theorem proven in |GN12] . 

Exponential convergence to equilibrium of the geometric fiber lay-down process 
towards a unique stationary state with explicitly computable rate of convergence. 



see Theorem 6.4 Here we use modern Hilbert space methods from the theory of 



hypocoercivity developed in [DMSIOJ . 
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2. Setup and Notations 

Before starting we introduce some notations, fix the language concerning the d-sphere 
and give a short introduction to the concept of Stratonovich stochastic differential equa- 
tions (SDEs) on manifolds. Until the end of this article we follow the notations and 
language introduced in the underlying section without further mention this again. 

C°° (X) denotes the set of all infinitely often differentiable functions / : X — )• M on some 
differentiable manifold X. The index c indicates compact support. V (or also denoted 
by V^ or Vjgd) always denotes the usual gradient operator in M'^, d G N, (with respect 
to the variable x) as column vector. | • | is the standard euclidean norm. The standard 
euclidean scalar product is simply denoted by • or also by (•, Oeuc- Superscript T denotes 
the transpose of some matrix. The expression smooth means that the underlying object 
is of class C°°. / is the identity matrix. Partial derivatives in M with respect to some 
variable x are denoted as usual by ^ or for short by dx- Convention: Any vector x G M 
is always understood as column vector. And the notation {x,y) for x £ M*^, y S M*^, is 
understood as column vector. 



Next let us fix some geometric language. Let d £ N. We consider the d-sphere given 
by 8''' = {v £ ]R'^+^ I |wp = 1}. The algebraic tangent space T^S'^ at the point v £ 8'^ can 
naturally be embedded into M^"''-'^. Under this identification any M-derivation D £ Ty8'^, 
V £ 8 , corresponds to some A £ M "'"^ with {A,v)^^^ = 0. We canonically identify D and 
A, in notation D = A. For / e C°°(S'^) it holds Df = A- V^d+if{v) where / is any 
smooth extension of / defined in an open neighbourhood in M'^^^ of v. So this justifies 
the notation Af = A ■ V^/(v) instead of Df. Moreover, for a given smooth vector field 
A on 8'^ defined by S'^ 9 u h^ A{v) £ Ty8'^, we write Af{v) = A{v)f, f £ C°°(S'^). The 
spherical gradient of such an / is denoted by grad§d/. One has 

grads,/(i;) = {I - v <S, v) V^,+if{v), f £ C°°(S'^), v £ 8''. 

Here x y := xy , x,y £ M.^^ and / is chosen as before. This definition is again 
independent of the local smooth extension / for /. In short notation we also write 
giadsdfiv) = {I-v(S)v) V^f{v) for / £ C°°(S^) and v £ 8'^. 

2.1. Stratonovich SDEs on manifolds. Now let X be a C°°-manifold, assumed to be 
second-countable and Hausdorff. X = X U {A} is the one-point compactification. Our 
main reference is jHsu02] . Consider also the excellent german book [HT94j . 

Solution concept. See |HT941 Def. 7.41] or |rW89l Ch. V]. Let (Q, J",P, {J"t}t>o) be a 
standard filtered probability space equipped with an r-dimensional standard {J^t}t>o- 
Brownian motion W = {Wt}t>o- Let Vq, Vi, . . . , Vr be smooth vector fields on X and let 
xo : ^ — )• X be J-Q-measurable. A solution X = {Xt}t>o of the Stratonovich stochastic 
differential equation 

r 

dXt = Vo {Xt) dt + Y, V, {Xt) o dWt^'^ (2.3) 

i=i 

with initial condition Xq = xq is any {J^i}(>o- adapted, continuous process on X having 
A as a trap such that the following is satisfied: Xq = xq P-a.s. and for every / £ C^(X), 

fiXt) - /(Xo) = / (Vo/)(X,) ds + J2 f (V,/)(X,) o dW,^^\ t > 0. 
Jo ^7o 
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We call e{X) := infj>o{^i = A} the explosion time or lifetime of X. Herein /(A) := 
for / G C^(X). Note that the definition of a solution X is always understood relative to 
(il, T, P, {Tt}t>o, W). In short form we write: X is a solution to SDE(Vo, • • • , Vr] W, xq). 
Finally, consider |Hsu02| Def. 1.2.3] for an equivalent definition involving test functions 
from C°°(X). 

State space transform. See | Hsu021 Prop. 1.2.4]. Let X and X be diffeomorphic with 
diffeomorphism ip : X ^ X. Then X = {Xt}t>o solves SDE Q iff X = {Xt}t>o, 



where Xt := (f{Xt), solves its associated SDE on X in which the Vj are replaced by their 
corresponding pushforward vector fields Vj on X. Clearly it holds {Vjf){p) = Vjf{p) where 
f o (p = f, p = ^{p) for / G C°°(X), p G X. The SDEs on X and X are called equivalent. 

Generator of a Stratonovich SDE. See |Hsu02| Sec. 1.3]. Any solution X = {Xt}t>o of the 
Stratonovich SDE ( |2.3| ) is an L-diffusion process generated by the second order Hormander 
type operator L : C°°(X) -^ C°^(X) given by L = Vo + lT^=iV]. Two such solutions 



X, Y (possibly defined on different probability spaces) of (2.3) with the same initial law 
weakly coincide, i.e., they induce the same law on the path space VF(X). Here W(X) is 
the set of all continuous paths u) : [0, oo) — t- X satisfying a;(s) = A for all s > t whenever 
uj{t) = A. W(X) is equipped with the canonical filtration. 

Brownian motion on Riemannian manifolds. See |Hsu02| Sec. 3.2]. In case (M,5() is a 
Riemannian manifold, any ^ AM-diffusion process is called a Brownian motion on M. Thus 
a Brownian motion can be generated by a Stratonovich SDE whose generator is equal to 
|Am. Here Am denotes the Laplace-Beltrami on (M,5f). 

3. The fiber lay-down geometry 

The three-dimensional fiber lay-down model is already presented in JKMW12] . Consider 
the latter for further interpretation. In this section we give a second view by view. As 
described in the introduction we now present a new differentialgeometric derivation and 
formulation of the model in each dimension motivated by the original ideas from |KMW12] . 
By the way, this shows up a fascinating interaction between applied and pure mathematics 
and illustrates the highly geometric nature of our model. We make use of the notations 
and the geometric language introduced in Section [2} For further background information 
in stochastic geometry, we again refer to |Hsu02] . 

3.1. The two-dimensional fiber lay-down SDE on M^. We start with the basic two- 
dimensional model from [GK MW07] which describes the lay-down of a single fiber as a 
curve ^t, t > 0, in the two-dimensional plane. The model is formulated in the introduction. 



see Equation (1.1). Therein, note that adWt = aodWt. Here o signifies the Stratonovich 
integral. 



Let X = {Xt}t>o be the strong solution of SDE (1.1 ) associated to some fixed standard 
one-dimensional Brownian motion W = {Wt}t>o- We use the notation Xt = (Ct,5t), 
t > 0. In Sectional we will see that a strong solution solution exists whenever $ G C°°(M^) 



although the drift vector of the SDE is not globally Lipschitz continuous in this case. 



3.2. The two-dimensional model on M? x T. The third component of the SDE (1.1) 
has the interpretation of being an angle. Henceforth, it is natural to formulate the previous 
SDE on M^ X T with T := ]R/27rZ. This can be done as follows. Consider the canonical 
projection P : M^ — ;► M^ x T mapping (^,a) to (^,a) where a := [a]. Then we introduce 
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the process X = {Xt}t>o defined by Xt := P(X(), t > 0. Herein X is the strong solution 
associated to (1.1) defined previously. We write Xt = (^t,at), at = [5*]. We claim that X 
solves the Stratonovich SDE 

dXt = AQ{Xt)di + Ai{Xt)odWt (3.4) 

on the manifold M^ x T with vector fields ^O) -^i given by 

A:=.-V,-V,*..-|-.^,:=.|-. (3.5) 

Indeed, let / G C;?^(M2 x T). In particular, / := / o P g C°°{W^). Furthermore, note 

vc/(p) = v5/(p), ||(p) = |^(p), p = (^,5) e^\p = nv). 

Thus the Stratonovich transformation rule (see e.g. |Hsu02| ) yields the claim since 
f[Xt) - /(Xo) = fiXt) - RXo) 

= J {T{as)-Vj{Xs)-yH^s)-T^{as)'^{Xs))ds + aJ ^{Z) o dWs 

= f {Aof){Xs)ds+ f {Aif){Xs)odWs. 
Jo Jo 

Summarizing, Equation (3.4) may reasonably be seen as a natural formulation of the 

two-dimensional fiber lay-down model on M^ x T. 

3.3. The equivalent tw^o-dimensional model on M? x S^. Since T and S^ are dif- 
feomorphic with a o u, v = T(a), SDE (3.4) can equivalently be formulated on the 
submanifold M? x S^ of M^. To obtain the equivalent SDE on M^ x S^ we only have to 
compute the pushforward vector fields Aq, Ai on E? x S-*^ associated to the vector fields 
^0 and ^1 from (3.5). But first some notation: Each x € M^ is written in the form 
X = {^,v) G M2 X m2 ^^^ ^± ._ (_^2,ui)'^. Now for / e C^iR"^ x T) let /G C'=^{R'^ x S^) 
be given as 

m ■=/{?), p = {C,a)eR^xT, p = {^,v) = {tT{a))eR^xS\ 
So if ^ denotes the pushforward vector field of ^ we get 

^m = ^fip) = ^mr{a)) = T^{a) ■ Vj{p) = v^ ■ Vj{p) 
oa oa oa 

where we have used that any smooth function on the submanifold M^ x S^ can locally be 

extended to some smooth function defined locally around (^, v) in M^. Thus the vector field 

T^ on M^ X T corresponds to the vector field v • Vj, living on M^ x S^. The latter is again 

canonically identified with {0,v-^), in notation v-^ • V^, = {0,v-^), v G S^. Moreover, note 

that for 1! G S-^, the orthogonal projection Ilgi [v] : M^ -^ T^S^, Ugi [v] (y) = (/ — v v)y, is 

equal to n§i[u] (y) = y ■ f"*" v-^ since T^S^ = spanju-*-}. Recall, v v = vv^. Altogether, 

we get 

A = . • V, - (/ - . v)vm • V. -(_(,_, ^' ,)v$(e)) ' 

Ai = av • Vt, = (T ( I 
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and the equivalent Stratonovich SDE on M? x S-"^ associated to (3.4) reads 

dXt = MXt) dt + MXt) o dWt (3.6) 

where W = {Wt}t>o is a standard one-dimensional Brownian motion. 

3.4. The d-dimensional model on M'^ x S"^"^. In the following let d € N with d > 2. 
Next, we translate the fiber lay-down model (3.6 1 in the most natural way to M x S . 
The resulting equation is then called the d- dimensional or geometric fiber lay- down model. 
But first note that the stochastic part of (3.6), i.e., the term ^i o dWt, models Brownian 
motion (with diffusion constant a) in the components representing S-^. This holds since 
the vector field V := v^ = v^ ■ V^2, v & S^, satisfies V^ = A§i. Now there is another 
possibility of writing A§i and generating a Brownian motion on S^. Therefore, set 

Vj := {I — V ^ v)ej = {I — V ^ v)ej ■ Vig2, v G S""", j = 1,2, 
and Cj being the j-th unit vector. Then one easily verifies J2j=i ^? = ^si- 

Having the latter identity in mind, we introduce the d-dimensional model, d > 2, on 
the manifold M*^ x S'^'^ as 



dXt = MoiXt) dt + Y^MjiXt) o dW^ 



U) 



(3.7) 



where W = {Wt}t>o is a standard d-dimensional Brownian motion and the vector fields 
J\fj are defined on M'^ x S'^^^ via 



AAo 



AA, 



V • Vg - (/ - w t;)V$(C) • V^ 
= a{I- 



I vjCi • V„ = a 







-{I-v0v)V^{C) 



j = l,...,d. 



(3.8) 



and any point from M x S is written in the form 



Here ej is the j-th unit vector in 
(^,f). $ G C°^{R'^) is only depending on ^. In case d = 2, SDE ( [3?7| ) then sHghtly differs 
from the original equation derived in (3.6). Nevertheless, their generators are the same. 
In particular, for d = 2, any solution to (3.6) coincides weakly with any solution to (3.7). 



Altogether, Equation (3.7) may reasonably be seen as the natural d-dimensional version 
of the fiber lay-down model. The physical relevant scenario are considering the cases d = 2 
and d = 3. This abstract SDE can now be embedded into M^ by extending the vector 
fields J\fj arbitrary to smooth vector fields on M^ . The extensions are still denoted by Mj. 
The solution to (3.7) is then obtained by solving the extended Stratonovich SDE in M^*^, 
see |Hsu02| Prop. 1.2.8]. By choosing the trivial extensions, we may consider the following 
Stratonovich SDE in M^'^ of the form 



d6 
dvt 



vtdt 



(3.9) 



-(/ -vt® vt) V$(Ci) di + (j{I -vt® vt) o dWt. 



Its solution stays on M'^ x S"^ ^ provided that the initial value lies on the manifold and 
gives a solution to (3.7), see |Hsu021 Prop. 1.2.8]. Here we define 



{I-vt(g) Vt) o dWt := ^(/ -vt(E) vt)ej o dWf 



U) 



j=i 
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3.5. The associated Kolmogorov operator. Let us calculate the Kolmogorov operator 
L : C°°(M"' X §"'-!) ^ C°°{W^ X S"'-!) associated to ([3?7j). As in the two-dimensional case 
it holds J2j=i-^f = ^^'^gd-i) see |Hsu02| Theo. 3.1.4]. Thus the associated generator 
L = Mo + l Tfj=i A/"/ is given by 

2 

L = vV^-{I-v(^ v)V^{i) • V^ + y Agd-1 (3.10) 

2 

= vV^- gradgd-iF + — A§d-i 

with y : M"' X S'^"^ — )■ M, V{S^,v) := V$(^) • v. Let us already remark that a stationary 
solution to the associated Fokker-Planck equation is explicitly known and is given up to 
normalization by e~^ ' , see Section [5| Convergence to equilibrium of the process X, 
solving the d-dimensional fiber lay-down SDE, towards its stationary state is studied in 
Sections [5] and El 



Remark 3.1. The generator L : C"^(]R^ x T) ^ C"^(]R^ x T) of (3.4) reads 



ft 2 a2 2 

L = r-V5-Vc^(e)-r^^ + y^ = r.V5-gradTrF+yAT. 



Now y : M^ X T -^ M is defined as V{^,a) = V$(C) • r. Thus we see that ( |3.10D is 
the natural generalization of the two-dimensional fiber lay-down generator from Mrx T. 
So, in order to derive our d-dimensional equations, we could alternatively translate the 
generators to higher dimensions and search afterwards for a corresponding modeling SDE. 



Following this approach we end up again with the SDE (3.7) 



3.6. Numerical simulations. For simulating our previously defined d-dimensional fiber 



lay-down SDE (3.7), we may of course simulate SDE (3.9) directly in the underlying 
euclidean space W^. Nevertheless, this requires a consistent numerical algorithm staying 
on the manifold M"* x E>'^~^. 

Instead, we take another approach and rewrite our fiber lay-down SDE in local coordinate 
form. Therefore, let d £ N, d > 2, and recall the following spherical coordinate system 
Td-iiO) = Td-i{6i, . . . , 6d-i) with 9i £ (0, 2tt) and 9j £ (0, vr) for j = 2, . . . , d - 1, see the 
Appendix. Furthermore, note that Td-i is 27r-periodic in 6i. Consider the following SDE 

dC = Td^i{e)dt (3.11) 

dOj = l-Gjie) v$(0 • nj{9) + y g]{9) {j - 1) cot(0j) J dt + agj{e) dw^^^ 

with j = 1, . . . ,d — 1. For abuse of notation, the time index t is omitted. Here Qj (with 
Qd-i := 1) and rij are given as 



and W denotes a standard {d— 1) -dimensional Brownian motion. Note that (3.11) writes 
as Stratonovich SDE in the same form. Let us justify this definition in the upcoming 
remark. 
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Figure 1. Two-dimensional case 



Remark 3.2. SDE (3.11) should be understood analogously as in Section 3.2 as some 
manifold-valued Stratonovich SDE with state space M x T x (0,7r) , or state space 
M^ X T in case d = 2, respectively. By using Formulas (7.27) and (7.28) from the Appen- 

he gene 
generator L computed on 



dix, observe that the generator corresponding to (|3.11|) coincides with the fiber lay-down 

^"^ X T X (0,7r)'^- 



Now assume that there exists a solution 
-Id-i) and having infinite lifetime. Then it is easy to 
T(i-i{Oi{t)^ . . . ,Od-i{t)))., is a L-diffusion. Hence X 



Y to (3.11) starting from (^,0i, . . . , 

see that X, defined by Xt := {i{t). 

coincides weakly with any solution of the (i-dimensional fiber lay-down SDE (3.7) starting 

from (^,rd_i(6'i,...,6'd_i)). 



Summarizing, the simulation of SDE (3.11) gives us a (local) L-diffusion, L being our 
d-dimensional fiber lay down generator. Thus it is reasonable to call (3.11) simply the 
d- dimensional fiber lay down SDE in local coordinate form. Note that in case d = 2, 
(3.11) reduces to (1.1) (or (3.4) respectively) and in case d = 3 we obtain back the three- 
dimensional fiber lay down model derived in |KMW12] . In Figure [T] and Figure [2l let us 
illustrate and compare the ^-trajectories in case d = 2, d = 3, for different values of a 
where a = 0, 0.1, 0.5, 4.0. $ is chosen as $ = |^p. 
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Figure 2. Three-dimensional case 



4. Global solutions to the stochastic equations 

Before studying the long time behaviour of solutions to our d-dimensional fiber lay- 
down equations, we shall discuss basic existence statements of the underlying stochastic 
equations itself. 



4.1. The two-dimensional model. First we analyze the basic two-dimensional model, 
i.e., the Ito SDE (lT^ in M^ <I>(0 = a^f + 6^|, o, 6 G M, a, 6 > 0, treats the physical 



relevant situation, see |KMW09] . It is easy to verify that in this case the drift vector of 
the underlying SDE is not globally Lipschitz continuous and one may ask under which 
conditions the stochastic equation admits a global solution. Therefore, we rewrite the 
equation in a bit more general form as 



d^i = T{at) dt 

dat = ^(6) • T^iat) dt + adWt. 



(4.12) 



The following proposition shows that even under the assumption of ^ being locally 



Lipschitz, explosion is not possible and SDE (4.12) admits a classical globally defined 
strong solution. 
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Proposition 4.1. Let ^ : M^ — > M^ 6e locally Lipschitz continuous. Given a filtered 
probability space (J7, J^, P, {J^t}(>o) together with an one- dimensional {Ft}t>o-Brownian 
motion W = {Wt]t>o- Then for each J^Q-measurable xq : ^ ^ M.^ there exists a unique 
{J^t}t>o-o-dapted continuous R.'^-valued process X = {Xt}t>o which satisfies Xq = xq F- 
a.s. and solves (4.12) in integral form. 

Proof. The classical Picard-Lindelof iteration scheme can be applied pathwise for each 
u; G J7. This is due to the special structure of our equation. We follow the strategy from 
p^S9n Ch. 5; Ex. 2.19]. So let t > 0, n G Nq = N U {0}, and define 



X 



(0) 



xo, X,("+^^ := xo + / 6(X(")) ds + AWt 
Jo 



where b{C,a) := (T(a), ^(0 • r^(a))^, A := (0,0, cj)^. Then X(") = {Xf>}t>o, n e No, 
is easily seen to be continuous and {J-j}t>o-adapted. 

Now fix T > and uj ^Q. Observe that the first two components of X^ (w) are bounded 
by N := \xo{ijj)\ + T for each n G No, t S [0,T]. Now N only depends on u and T. 

Thus we may restrict b to [—NjN]"^ x M in the definition of X^ (co) for all n E No, 
t E [0,T]. Furthermore, observe that b is Lipschitz continuous on each K x M., K C E? 
compact. In particular, there exists some constant L, depending on to and T, such that 



\b{x) -b{y)\ < L\x-y\ for all x,y E [-iV,iV]- 



I. So altogether 



6(xi")M)-6(x("-^)H) <L x[-\u:) 



Xt'\u] 



(4.13) 



holds for all n E N, t ^ [0,^]. Now we may continue as in the classical case. Indeed, for 
those n and t we define 



Dt\^) 



max 

0<s<i 



Xi")(..)-X("-1)(..) 



i(")/ 



^{n-l). 



Then dils] ) yields DI'''{uj) < L J^ D'-J' ''{u:)ds' for n > 2 and t E [0,T]. 

d[^\uj) is bounded for t E [0,r] by K := T\b{xo{uj))\ + a maxo<s<T \Ws{oj) 
depends only on cj and T. Thus inductively we get 

(Li)"-i 



Moreover, 
which also 



Dt\u:)<K 



(n-1)! 
Now let n > ?7i, t E [0,T]. The latter estimate implies 



riEN, tE [0,r]. 



max 

Q<s<t 



xW(a;)-X(")( 



OJ, 



ri-l 

< > max 
^^ 0<s<t 
]=m 



Xp+i)(a;)-Xi^-)(a;) < K Y, 



(Lt) 



Hence X^'^'{uj), n E N, converges uniformly on [0, T] to some continuous function Xt{uj) 



where < t < T. Further note that Xt{u}) is contained in [—N, N] ■ 



Rfor all t E [0,r], 
since Xf"''{co) satisfies this property independent of n E N. Hence for t E [0, T] we conclude 



6(Xi")(a;))ds- / 6(X,(w))ds 



<tL max 
0<s<t 



Xi-\u) - X.ico] 



0. 



This shows that Xt{oj), < t < T, solves (4.12) in integral form. Now the previous 
construction holds for all T > and each uj £ 0,, hence the map 1 1— >■ Xt{uj), t > 0, is well- 
defined and continuous. Clearly, Xt is J^(-measurable, t > 0. Concerning the uniqueness 
statement, see e.g. |KS91[ Ch. 5, Theo. 2.5]. D 
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Now we are able to solve easily the two-dimensional fiber lay-down SDE (3.4) with 
state space R^ x T without using any abstract arguments. For simplicity we shall assume 
$ G C°°{M?) to stay consistent with the (smooth) manifold language. Concerning all the 
following uniqueness statements, we refer to |Hsu02| Prop. 1.2.9]. 

Corollary 4.2. Let ^ £ ^^(M?). Then there exists a unique solution X = {Xt}t>Q with 
infinite lifetime of the two-dimensional fiber lay-down SDE(Ao,Ai;W,xo), see Equation 
(|3.4p, with state space M^ x T. 



Proof. Consider the canonical projection P : M^ — ^ M^ x T. Now we may choose some 
xo : il — ^ M^ being J^o-^ieasurable that satisfies xq = P o xq. Such an xq exists since the 
map (f : [0,2tt) — )■ T, x i— ^ [x], is a Borel isomorphism, see |Par67| Cor. 1.3.3]. Proposition 
4.1 (with ^ = — V$) is applicable and yields the existence of a strong solution X to SDE 



(4.12) with initial condition xq relative to {Q,J^,¥,{J't}t>o,W). The process X, defined 



by Xj = P o Xt, t > 0, satisfies Xq = xq P-a.s. and has infinite lifetime. The fact that X 



solves SDE(^0)-4i; W,xq) has already been verified in Section 3.2 D 



4.2. The d-dimensional fiber lay-down model. The same statement as before remains 
true for our general d-dimensional equation. 

Proposition 4.3. Let ^ G C°^{M. ), d > 2. Then there exists a unique solution X with 
infinite lifetime to the d- dimensional fiber lay-down SDE{J\fo, . . . ,Mn', W,xo), see Equation 



(3.7), with state space M'* x S 



d s, sjrf— 1 



Proof. By |Hsu02| Theo. 1.2.9] we know that there exists a unique continuous solution 
X to SDE(A/'o, . . . ,J\fn', W,xo) with state space R^ x S'^~^ up to its lifetime e. It remains 
to check that e{X) = oo holds P-almost surely. Write Xt = {Ct,vt), t < e{X). For each 
/ G C^{R'^) we have 

/(6) - mo) = fvs- ^di^s) ds,0<t< e{X). 
Jo 

In particular, this identity holds for all / G C^{M. ) which satisfy /(^) = ^j for j = 1, . . . ,d, 
inside some balls -Br(O) of radius r > large enough. Then we easily conclude 

^t-^o= f Vsds, 0<t<e{X). (4.14) 

Jo 

Consequently, |Ct — ^o| < t for < t < e{X). Now assume that e{X) < cx) holds on some 

A £ T with F{A) > 0. But then |X(|]jj2n — ^ cxd and hence |^t| — )■ oo as 1 1 ^i^) on A, see 

e.g. JHsu021 Prop. 1.2.6]. But this contradicts the previous inequality. D 

5. Convergence to equilibrium: Mixing properties 

As mentioned in the introduction, we start in this section with the investigating of the 
long time behaviour of our d-dimensional fiber lay-down process by proving strong mixing 
properties. We make use of a new version of Doob's theorem derived recently in |GN12] . 
But first we introduce some notations that are used until the end of this article. 

5.1. Motivation and setup. In the following let d G N, d > 2. We focus attention on 
the d-dimensional fiber lay-down SDE with state space M := M^ x S'^^'^ 

d 

dXt = Mo{Xt) dt + J2^jiXt) o dWt^^^ (5.15) 
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where VF is a standard d-dimensional Brownian motion. For computational convenience, 
we replace <I> by 3^*^ in SDE (5.15). So Ao is redefined as 



Mo{^,v) 



-:j{I - V (S> v)V<^{^) 



{^,v)e 



l'^ X S'^- 



p2d 



The vector fields Mj are defined in (3.8). Convention: The first d-components of M'" are 
abbreviated by S, whereas the last d- variables are denoted by v. We consider potentials 
$ : M*^ — 7- M such that e~* G L^iJR.'^, d^). Thus without loss of generality we assume in the 
following 

'de = l. 

,13{M)) defined by 

' 5. 



We introduce the probability measure ^ on (} 



/^ 



^dS, ^ v, v :- 



-iF+^A, 



IS 



vo^S-^-i) 

Here S denotes the surface measure of S , vol(S ) the surface area of S and ;S(M) 
the Borel-sigma-algebra on M. The Kolmogorov operator L = w-Vg— gradgd 
see (3.10), associated to SDE (3.7) is also sometimes denoted by L^ where V : 
now given as V{^, v) = ^^^-j- V<^(^) • v. 

Moreover, we denote by X^ = {Xf}t>o the solution to SDE (5.15) with initial value 
X G M. Formally, its corresponding probability density ft (with respect to d^ (8) 5 and 
Dirac initial state x) satisfies the Fokker-Planck equation dtf = L f with L being the 
Fokker-Planck operator associated to the d-dimensional fiber lay-down model. The latter 
is given as the (algebraic) adjoint of L in L^(M,d^ S). Thus, using Lemma 7.2 from 
the Appendix below, we get 

,2 



L 



FP 



a 



Vt + gradgd-iF - V$(0 • t; + — ^ 



So, up to normalization, a stationary solution to the Fokker-Planck equation is now given 
by 

-HO 



Hence we expect that ft converges towards 



volC S''-M -^ as t — !► cxD. In other words, Xf should 
be distributed accordingly to ^ for large values of i > 0. Now formally, n(t, x) = E[uo(Xf )] 
solves the Kolmogorov PDF dtu = Lu with initial state u{t = 0) = uq. Hence altogether 
the subsequent analysis in the previous and the upcoming section has the interpretation 
of proving convergence of E[uo(-'^f )] towards E[no(X^)] := J^uodfi. We now start with 
the mentioned strong mixing properties. 

5.2. Mixing properties. First, we need the following lemma. Therefore, for given vector 
fields Vi , . . . , Vr , on some manifold X we denote the least M- vector space including all Vj , 
j = 1, . . . , r, which is closed under the Lie-bracket operation by L(Vi, . . . , Vr). 

Lemma 5.1. For the vector fields Mj, j = 0,1, . . . ,d, from (3.8) it holds 

d 



dim L{J\fi, . . . ,^fd,^fo + 



dt^ 



2d at each point of (0, 00) x 



where we assume a > 0. 
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Proof. We may set a = 1 and choose p = {t, ^, v) £ (0, oo) x M arbitrary. First note that 
span{v} e span{(I - -u v)ej \ j = 1, . . . ,d} = span{v} T^S'^^^ = M'^. 



Observe that 



AAj,AAo + 



A 



J •- 



Mj,Mo + 



for J = 1 , . . . , d, is of the form 
d 



dt 



for some smooth function f-^ 
Appendix, we obtain 

d 



U) 



M'^. Hence, by the proof of Lemma 



7.1 



from the 



^:=E 



i=i 



M, 



AAj,AAo + 



dt 



-{d-l)vVt: + f2{^,v)-V, 



Here /2 : 



.'^ is again smooth. So finally, under 
5 



AAo + 



dt 



(p), AAi(p),..., AAd(p), ^i(p),..., Ad{p), A{p) 



we may always choose 2d-linear independent vectors. The claim follows since the manifold 
(0, oo) X M has dimension 2d. D 



We need one more lemma. X'^ denotes the process solving (5.151 with state space 
M defined on some underlying probability space (Jl, J^, P, {J-i}i>o) satisfying the usual 
conditions. Br{z) denotes the open ball with radius r > centered at z £ M . 

Lemma 5.2. Let $ G C'^iR'^) and let f G C^(M). The function u{t,x) := E[/(Xf)], 
t > 0, X G M, is continuously differentiable in t, twice continuously difjerentiable in x and 
satisfies 

d 

—u{t,x) = Lu{t,x) 

where L : C°°(M) -^ C°°(M), L = Mq + I E?=i AA/. 

Proof. This is a local statement, so let x = {^^,v^) £ 

fixed. Recall that the solution X^ = {Xf}t>o, y G IV 

i.e., SDE(A/'o, . . . , A/^; W, y), can be obtained as the solution to the embedded Stratonovich 

SDE in M. in which the vector fields J\fj are extended to M? in the obvious way. Define 

ri = r, r2 = 2T, and choose some Xi;X2 G C^{M.'^) with xi = 1 on Br2{^^), X2 = 1 on 

-62(0), < xi, X2 < 1 and define Afj on M^o! as 

!)2d 



II and let T > be arbitrary but 
to the abstract equation (5.151, 



^M,v)= xi{0X2iv) MjiCv), {C,v) G 



j = 0,...,d. 



Now let X^ be the solution to the Stratonovich SDE 

d 



dXt = AAo dt + ^ Mj o dWt 



U) 



(5.16) 



m M""- with initial condition Xq = y £ M. All Mj are still tangential to M. Thus by 
|Hsu02| Prop. 1.2.8] we conclude that X^ stays on M and, as in the proof of Proposition 



4.3 



Xy has infinite lifetime. Hence for each y = {^^ jV^) G M with |^^ — ^^'| < ri it follows 
that Xy and Xy both solve the Stratonovich SDE ( |5.16[ ) in R^d ^p to time T. Clearly, the 
latter SDE can be written in Ito-form with global Lipschitz coefficients. Consequently, we 
have X^ = X^ P-a.s. for all < t < T and all y as specified above. 
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Now choose any extension of / to some / E C^{M. ). And since the equivalent Ito-form 



of (5.16) has C°°-coefficients having compact support, we get that u{t,z) := E[/(Xf)], 
i > 0, z G M. , is continuously differentiable in t, twice continuously differentiable in z 
and satisfies 

d ^ 

-Q^u{t,z) = Lu{t,z) 

where L : C°°(M2'i) ^ C°^(m2"!), L = Mo + E?=i AA/, see |DPRlll Prop. 2.7] and JBF6T] . 
Hence the claim follows since u{t,y) = u{t,y) for < t < T, y G i?ri(^^) x S*^"^, as well 
as since Lh{y) = Lh{y) holds for each h E C^(M) and all such y. D 

We remark that under the assumption of the previous Lemma, u{t, x) is even infinitely 
often differentiable in x, see |BF61| . Nevertheless, we do not really need this stronger 
conclusion. Now we end up with our desired theorem. i?(M) (or C(M) respectively) 
denotes the set of all Borel-measurable (or continuous respectively) real-valued functions 
on M. The index b denotes all bounded functions of the underlying set of functions. 

Theorem 5.3. Let d e N, d > 2, and let $ G C°^{W^) with J^a e''^ d^ = 1. Assume 
that a > 0. Then the d- dimensional fiber lay down process is strongly mixing, i.e., for all 
f £ Bb(M.) we have 

limE[/(Xf)]= / /d/i 

uniformly m x G M on compact subsets o/M. 

Proof. We apply a version of Doob's theorem, see |GN12| Theo. 4.6] . So we define oper- 
ators Tt, t > 0, Tt : Bb{M) -^ Bb{M) as Ttf{x) = E[/(Xf )] for / G Bb{M) and x G M. 
These operators are indeed well-defined. Therefore, let X be the process defined on the 
path space W{M) by Xt(w) = u{t), t > 0, w G W{M). Then it holds E[/(Xf )] = E^[f{Xt)] 
where Kx denotes expectation with respect to P^ = Fo[X^)~^. Thus for A G ;S(M) and all 
t > we have Ttl^ix) = ¥x{Xt G A} which is i3(M)-measurable in the variable x G M, see 
|Hsun2l Sec. 1.3]. Hence Ttf G ^^(M) for ah / G ^^(M). Furthermore, each Tt is a Hnear, 
positive operator satisfying ||Ti/||oo < ||/||oo for all / G Bb{M.). Here || • ||oo is the usual 
supremums norm. Moreover, for A G B(M.) and x G M we have TtTslA{x) = Tt^glj^^x). 
This follows by the strong Markov property of the L-diffusion X^, see |Hsu021 Theo. 1.3.7]. 
By applying monotone convergence we obtain TfTgf = Tt-\-sf for each / G Bb(M.), i.e., 
{Tt)^^Q is a semigroup. 

Nextlet / G C^{m). The function u{t,x) := E[/(Xf)], t > 0, x G M, is continuously 
differentiable in t, twice continuously differentiable in x and satisfies 

d 

—u{t,x) = Lu{t,x), 



see Lemma 5.2 Now choose r > such that the support of / is contained in Br{0) x S 



d-l 



By using Identity (4.14) observe that for fixed T > 0, the function n(t, •) vanishes outside 



K := Br+T{0) X S^-^ for t G [0,T]. Moreover, note that ^'u(t,x) = E[L/(Xf)] holds by 
Ito's formula for all t > 0. Thus 



-u{s,x) 



< sup \Lf{y)\ < oo, s£[0,T], x£ K. 
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Note that ^u{-, •) is 'B[0,t] ;B(M)-measurable for all t > 0. Now integrate the identity 

^(*; ■) ~ /(■) = /o ^^(•^j ■) 'i'5' ^ ^ [0' ^]' with respect to //. By the previous we may apply 
Fubini's theorem and get 

f Ttfdfi- f fdfi= f f Lu{s,x)dfi{x)ds = 0, te[0,T], 
JM Jm Jo Jm 

where the last equality follows since J^^ Lhdfi = holds for all h € Cf (M), see Proposition 



7.2 from the Appendix below. Hence we have J^Ttf d^ = J^fdn for all / G C^(M). 
Thus by the functional monotone class argument, see |BG68| Ch. 0, Theo. 2.3], we easily 
conclude Jf^ Ttf d^ = J^^f d^ for all / G Bb{M.), i.e., fi is an invariant measure for (Ti)t>o. 
Furthermore, Lemma |5.1| just says that condition (E) from |IK74j is satisfied. Hence 
by tIK74i Theo. 3] we obtain that there exists pt{x,y) £ C°°{{0,oo) x M x M) with 
P{Xf G dy} = pt{x,y)dy, t > 0. And since Xf has infinite lifetime, i.e., TjIm = 1m 
holds for all t >0, the existence of such a function p then implies by [IK 74' Lem. 5.1] that 
{Tt)t>o is already strongly Feller. Here strongly Feller means TtBh{M.) C Cfe(M), t > 0. 
Next note that condition (ii) of Propositon 4.3 from |GN12j is obviously satisfied since n 
is stricly positive and 1m plays the role of the element e in |GN12| Prop. 4.3]. Together 
with our previous conditions this already implies the weak irreducibiliy assumption from 
|GN12| . 

It is left to verify the strong continuity assumption in L^(Wl,fi), i.e., condition (4.1) in 
|GN121 Theo. 4.6]. To verify the latter, it suffices to show that limj^o Ttf{x) = f{x) holds 
for all / G Cb(M.) and ah x G M, see |GN12| Rem. 4.7]. But this is obvious, use Lebesgue's 
dominated convergence. 

Finally, observe that also condition (iv) of |GN12| Theo. 3.4] clearly holds, hence the claim 
follows by the generalized version of Doob's theorem, see jGN12| Theo. 4.6]. D 

Remark 5.4. We remark that the original version of Doob's theorem, see |DPZ96j is not 
applicable in case of the fiber lay-down process. This is because the strong irreducibility 
assumption, i.e., the to-irreducibility property, is not satisfied in our case. The latter 
means that for some to > it holds Tt^lrix) > for all x G M and each non-empty open 
r C M. But we have already seen in the previous proof that such a property can't be 
satisfied simply due to the fact that the velocity vector of our fiber lay-down process lives 
on the sphere. 

6. Convergence to equilibrium: Hypocoercivity 

Now we switch to functional analysis. This section is independent from the stochastic 
ergodic results and statements obtained in the sections before. Now we apply the fasci- 
nating hypocoercivity theory derived in |DMS10) and generalize |DKMS12] to our higher- 



dimensional setting. For d = 2, Theorem 6.4 below contains the result from |DKMSl"2] . 



In contrast to |DKMS12| we stay in the Kolmogorov picture. In this way we obtain a nice 
comparision with |GK08| . However, due to the underlying L^-framework, the Kolmogorov 



setting is equivalent to the Fokker-Planck setting, see Remark 6.2, Hence we may switch 
between both. 

Remark 6.1. We remark that we do not specify any domains of all considered operators 
in this section. So our subsequent analysis in this section is done at first algebraically. 

6.1. Hypocoercivity. We closely follow [DMSIO', Sec. 1.3] and |DKMS12j . Remember 
the general setup and notations introduced in Section |5.1[ A convenient choice for the 
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underlying Hilbert space is 



L' 



>/^)> (5,Ml2(/.) 



M 



ghdfi, 



5.1 



For further motivation on this choice, 



with M = M"' X S^"^ d > 2, and /i as in Section 
the interested reader may also consult e.g. [GKQ8) or [CG08| . The corresponding norm 
on L'^{ij.) is denoted by || • ||i:,2(^). The Kolmogorov operator L^ = L associated to our 
d-dimensional fiber lay-down SDE is written in the form L = S — A where 



S 



1 



a^A, 



-v-V^+ grad§d-iy, V{C, v) 



1 
d-1 



v-^iO-v- 



Now Lemma 7.2 from the Appendix and the choice of jj implies that A is antisymmetric 
on L^(/i), S is symmetric and negative semi-definite on L?'{p) and we have 



/ Ludfj, 

Jm 



0. 



Finally, in the following, the Kolmogorov PDE 

dtu = Lu = Su 



Au 



(6.17) 



(6.18) 



is considered as an abstract Cauchy problem in L'^{ij,). Its solution subject to the initial 
condition uq £ L'^ifJ-) is denoted by u{t), t >0. Motivated by Section 5.1 we have to study 
exponential convergence of u{t) towards (uq, 1)2,2(^1 = J^uodfi in L^Jp) as i — )■ 00. We 



just remark that u{t), t > 0, when starting with (6.18), can really be given a rigouros 
meaning as some E^[«o(^t)]- This fact is justified by the theory of (generalized) Dirichlet 
forms, see e.g. |Tru05| . 

Remark 6.2. 

(i) Let us explain the equivalence between the Kolmogorov and the Fokker-Planck L^- 
setting. Following |DKMS12] . the underlying Hilbert space in the Fokker-Planck 

Now consider the Hilbert 



picture is just L?'{p) 
space isomorphism 



L\M,]1) with // 



'di 



V. 



T : l2(/x) -^ L^{Ji), h ^ Th, Th{i,v) 



-#(0 



= L^P. Hence n(t), t > 0, 



One readily checks that (algebraically) it holds TL T~ 
solves the abstract (Kolmogorov) Cauchy problem dtu = L u in Lp'{n) with initial 
condition u{t = 0) = wq if and only if /(t), t > 0, solves the abstract (Fokker- 
Planck) Cauchy problem dtf = L f in L?'{]1) with initial condition f{t = 0) = /q. 
Here f{t) := Tu{t), t >0. Moreover, note that 



u{t) -(no,l)i2(^) 



i.e., we have \\u{t) — Jj^ uq d^ \ 






f{t)-F{fo,F)^, 



M 



lHp) ' 



\\f{t) 



F /jyj /o d^dz^ 1 1 l2 (-r^ • This shows that 



Hm)- 



we may equivalently switch to the Fokker-Planck L^-setting in Theorem 6.4 below, 
(ii) By integrating the identity u{t) — uq = Jq Lu{s) ds with respect to ^ and using 



(6.17) we obtain 



{u{t),'^)LHu)= / u{t)d^i= / nod/i = (no,l)i2( ) 



(6.19) 



This has the interpretation that /x is an invariant measure 

(/0,-P)l,2(~) 



On the other hand 
that is, /m /(t) d^dz^ = /„ /o d^di/. So 



(I6.19D implies (/(t),F),,,~, 

in the Fokker-Planck L^-setting, this just means that mass stays conserved 
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Next we fix the conditions on the potential ^ analogously to [DKMS12] . Let d G N, 
d > 2. We denote by H'^{e~^d^), /c G N, the space of all /c-times weakly difFerentiable 
functions /i : M — )• M whose derivatives up to order k (including the function itself) are 
elements of L^(e~ d^). Further, V? or V^ stands for the Hessian matrix in M . 



2,00, 



(HI) Regularity: $ G W^^^iR"^). 

(H2) Normalization: J^d e~*d^ = 1. 

(H3) Spectral gap condition: There exists a positive constant A < 00 such that 



/ |Vu|2e-*de> A / 



uh~^d^ 



for ah u e Fi(e-*dO with /^d we"* d^ = 0. 
(H4) Pointwise condition: There exists a constant C < 00 such that 

V2$(e) <C{1 + |V$(OI) for ah C G M"^. 



Remark 6.3. Let d e N, d > 2. In analogy to |DKMS12] we need elliptic regularity 
estimates from |DMS10( Sec. 2] . This L^ — >■ H^ regularity result requires the full strength 
of all conditions on $. So assume (H1)-(H4) and let / G L^(e~*d^). Assume that 



u £ L^(e-*dO with /ig 



ue 



■de 



solves the elliptic equation 

--Au + -V$ ■\7u + u = f. 
d d 

Then jDMS10( Prop. 5] together with |DMS10l Lem. 8] implies the estimates 



|V «llL2(e-*dC) ^ C*! \\J\\L-i{e-fdii) 



l^l|Vn|||^2(,-*dO<C: 



L2(e-*dO 



with W 



1 + ^ |V$| . Here Ci, C2 < 00 are constants independent of / and u. More 

0^1 f- 



precisely, the elliptic equation from |DMS10| Sec. 2] is of the form ojq u—V-{ijj\ Vu) 



Note that the latter reduces for Wq 



e and oj( 



^e * to the one from above. Now 



in or der to ap ply |DMS10l Prop. 5, Lem. 8], Co nditions (2. 1), (2 .4), (2.5) , (2.6) and (2.7) 
from |DMS10] must hold true. Indeed, consider |DKMS12] and |DMS101 Sec. 3] for their 
verification. 

After this preparation, we end up with the final theorem. 

Theorem 6.4. Let d £ N, d > 2, and a > 0. Assume conditions (H1)-(H4). Then, for 
every r] > 0, the solution u{t), t >0, to the abstract Cauchy problem (6.18) in L'^{^) with 
initial condition u{t = 0) = uq G L'^ifJ-) satisfies 



u{t) -(U0,l)i2(^) 

Herein A is given by 

A 



LHt^) 



rj 



< 



(1 + 77) Uo-{uo,l)LHf,) 



Kia^ 



L2{m) 



and the constants Kj < 00 



j = 1,2,3, are only depending on the potential $. 

Proof. Step 1: In analogy to |DKMSl"2] we start by realizing the program from |DMS10( 
Sec. 1.3]. Therefore, we choose the desired Hilbert space Ti. as 



n 



g£L\M,^) (5,1)l2(;.) 



M 



gd/x = 
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Next, we define the deviation g{t) = u{t) — (uq, ^)L^(a) ; * ^ 0- By (6.19) it follows that 
git), t > 0, solves the abstract Cauchy problem (6.18) in the Hilbert space 7i subject to 
the initial condition g{t = 0) = uq — (uq, 1)^2/^) G Ti. Furthermore, denote the orthogonal 
projection to M{S) (the null space of S) by 



n^ = Pg-= / gdu. 

J§d-1 



In particular, J^d PgC * d^ = 0. Clearly, it holds AUg = —v ■ V^Pg and hereby n^II = 0. 
By |DMS10( Lem. 1], the latter identity implies 

\\Bg\\m^)<\\\{I-Il)g\\L2(^^^ (6.20) 

where B is given as i? = (/ + (An)*yin)~^(^n)* and I denotes the identity operator. 
For completeness, let us recapitulate the general strategy from |DMS10J and |DKMSl"2] . 
First, the modified entropy functional is defined as 

He[g] ■■= ^\\g\\h(^f,) + e (55,5)l2(^), eg (0,1). 



Then Hir[g] is equivalent to ||5|||2(„-), more precisely, (6.20) yields 



(m)' 
g\\h^^) < Hs[g] < ^—^\\g\\h(^^y (6.21) 



1 — e,, 1,9 ^T r 1 1 + e,, 1,9 

^-||5lli2(^)<i:^.[5]<^ ""'" 



Moreover, the evolution of Hi;[g{t)] is given by ^ H^[g{t)] = —Di;[g{t)] where D^[g] is the 
entropy dissipation functional 

DM = - {Sg, 5),.2(^) + e {BAUg, g)^,^^^ + e {BA{I - U)g, g)^,^^^ (6.22) 

- e {ABg, g)^2^^) - e {BSg, g)^2^^^ . 

The main step involves showing coercivity of D^] for some appropriate e G (0, 1). Indeed, 
assume this to be true. Then there exists k > such that D^] > k||(7|P and hence by 
(6.21) we have ^ Hi,[g{t)] < — j^i?£[g'(t)]. By using Gronwall's lemma and again (6.21), 
exponential convergence of g{t) towards in L^{p) as t — ^ cxd follows. The desired rate of 
convergence is then finally be determined in step 4. 
But first of all, we have to show coercivity Di;[-] for some e G (0, 1) small enough. 

Step 2: To do so, we start by verifying the microscopic and macroscopic coercivity as- 
sumptions from |iDMS10| Sec. 1.3]. The first one is satisfied due to the Poincare inequality 
on S"'-! 

- — - / (gradgd-i/i,gradsd-i/i)^55d-i dv > I h^ dv - \ I hdv] , 
d-l Jsd-i ^ Jgd-i \Jsd-i ) 

see |Bec891 Theo. 2]. This yields 

- (^5, 5)l2(;.) > \cy\d - 1) ||(/ - n)ff||i.(^) . 

Here (•, ■~)j,^d-i is the usual Riemannian scalar product on the tangent bundle of the 
sphere. Using Lemma 7.3 from the Appendix below, which is just a simple application of 
the Gaussian integral formula, and (113) we conclude 

Pn5||i.(^) = / / {v V^pgf e-" di^iv) de 



20 



M. GROTHAUS, A. KLAR, J. MARINGER, AND P. STILGENBAUER 



The latter inequality is the desired macroscopic coercivity property. Henceforth, see 
[DMSlOl Sec. 1.3], we have 

A 



{BAIig,g) 



L^(p) 



> 



d + A 



l|n<7ll 



L^ip) ■ 



So the sum of the first two terms in (6.22) is coercive. Now as described in |DMS10] and 
[DKMS12] , coercivity of D^ for e small enough follows if we can show that the operators 
BA{I — n), AB and BS are bounded and satisfy certain estimates. This is done next. 

Step 3: Again we proceed as in |DKMS12] . First of all, boundedness of AB is auto- 
matically satisfied due to our previous conditions, see |DMS10l Lem. 1]. More precisely, 
we even have 

P^5llL2(^)<IIU-n)5llL2(^). 

Furthermore, by using Formula ( |7.26[ ) from the Appendix, we conclude 

{AU)*g = -UAg = V^- f vg<lu{v)-j giad^d-iV{v)dv{v) 



(^^-^**)-i 



V gdv{v). 
This together with Equation (|7.25[) yields {AlifSg = -{d - l)^{AIifg. Consequently, 



BS 



-{d 



1)^B and therefore 



\\BSg\\ 



LHt^) 



<{d-l 



(7 



{I-m\ 



L2(/.) • 



Finally, boundedness of BA is equivalent to boundedness of (BA)*. To verify the latter, 
we make use of the elliptic regularity result mentioned before. First we calculate 



^iph v)^^^ - j—^ {{I-v® v) V^^, ViPh), 



(6.23) 



[vylph- 



d 



^ (Vc^, V^p,)^^^ + ^ {v, V^c^)^^^ {v, V^ph) 



Thus by the Gaussian integral formula, see Lemma 7.3, we get liA^Iih = ^A^p^ — ^V^$ 



V^Ph- For g £ L'^(M,p), let h = {I + {All)* {AIl))~^ The previous calculation yields 

Pg = nh- liA^Iih = -^^iPh + ^V^$ • Vt-ph + Ph- (6.24) 



Note that {BA)*g = Alih. Using (6.23) and applying the elliptic regularity result from 



|DMS10) (see Remark [63J) to Equation ( 

||(i?^)*5||L2(^)<|||V|/5,|||i2 

< C\\pg\\L\. 



6.24P , we get 

1 
+ 



(e-*dO ■ ^_i 



|V5^||V5P/,|||i2(,-*d^) 



'*dO < K\\9\\L'^{^l) 

for some constant K < oo independent oi g. Thus also ||i?^(7||2,2(„) < i^Hfi'llLa/^) and 

\\BAiI -U)g\\L2^^^ < K\\il -U)g\\L2^^y 

Step 4-' By jDMS10( Sec. 1.3] and the previous steps we infer coercivity of D^ for some 
e £ (0, 1) small enough. As explained at the end of step 1, this yields already exponential 
convergence. It remains to verify the claimed rate of convergence. To do so, just copy the 
calculation from subsection 3.4 in |DKMS12l Theo. 1] via replacing D, D := ^, through 
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^2 



D := (d — 1)^ and (2 + A) through d + A in the latter. Furthermore, note that in 
subsection 3.4 from |DKMSl"2l Theo. 1] the expression 1 + 6*2-0^ can be replaced by some 
term of the form 1 + C^D + C^D^ with C3, C4 being some finite constants not depending 
on (T. Having this in mind, the desired rate of convergence of the theorem is shown. D 

Remark 6.5. A solution to the (Kolmogorov of Fokker-Planck) Cauchy problem in case 
(J) £ C°°{M.) can easily be constructed using the well-known Heffer-Nier construction 
scheme, see |HN05| Prop. 5.5]. This result can then easily be generalized to the case of 
$ being locally Lipschitz continuous (and bounded from below) as in |CG10] . Details 
on this and further analytic approaches for proving convergence to equilibrium of the d- 
dimensional fiber lay-down process are discussed in a forthcoming publication of the first 
named and the last named author of the underlying article. 

7. Appendix 

For completeness, in this section we just prove some specific statements needed for our 
analysis in this paper. We start with the first one. 

Lemma 7.1. Let d E N. /§d denotes the function v ^-^ v, v £ S , where Agdigd is 
understood componentwise. Then it holds 

Agd/§d = -dl^d, (7.25) 

Proof. We make use of a specific representation formula for the Laplace Beltrami. So let 
Vj be the vector field on S'^ given by u 1— )■ Vj{v), with Vj{v) G T^S'^ the j-th row of the 
matrix I - v v, v e E,'^. It holds E?=! V| = Agd, see |Hsun2[ Theo. 3.1.4.]. Choose 
1 < J < (^ + 1- Now easily Vjlgd{v) = Vj{v). Thus the i-th component of Vjlgd{v) is just 
Pij{v) := 5ij — ViVj. We have 

^ d 



dVn 




-Vi + 2viv'^ if i / J 



if ^ = J 
and therefore 

d+l d+1 

^ Vj{Pij)iv) = J2 (-"i + ^^'i^'f) - 2i;i(l - vf) = -dv^. 

Hence we get Ag^Igd = Y.'j=i '^'j h^ = -dl^d. D 

The following proposition is used for computing the Fokker-Planck operator correspond- 
ing to our higher-dimensional fiber lay-down process. S denotes the Riemannian measure 
on (S'^,B(S'^)) and recall that C'^{E'^) is dense in L^(S'^,5). Here C^ means twice contin- 
uously differentiable. 

Proposition 7.2. Let d G N. V: S'^ -^ M is defined by V{v) = {z,v)^^^, where z £ M°'+^ 
is fixed. Define the vector field A by A = gradgdV and consider the operator L = A with 
domain C^(S'^). Then its adjoint operator L* in L^(S'^,5) is given on C^(S'^) as 

L* = -A + d{z,v)^^^ 



22 M. GROTHAUS, A. KLAR, J. MARINGER, AND P. STILGENBAUER 

Proof. Let f,g e C^S'^). Clearly, J^,A{fg)dS = J^aA{f)gdS + J^,fA{g)dS. So it 
suffices to show 

/ Aih)dS = df {z,v)^^hdSiv), h£C^{§'^). (7.26) 

So let h G C^(S'^). One has A{h) = gradgd{V){h) = {gTad^dV,gTad^dh)j,^d where (-JOTS'* 
denotes the canonically Riemannian scalar product on the tangent bundle TS of S . 
Hence we get 

/ A{h)dS= / (gradgdF, grad§d/i)^(jd dcS = — / {AgdV)hdS. 

Jgd J§d Jsd 

Finally, by Lemma 7.1 we have AgdV = —d(z,Igd)^^^. Thus (7.26) is shown and the claim 
follows. n 

Moreover, we need one more lemma, which is just a simple application of the Gaussian 
integral formula, u denotes the normalized Riemannian measure of S , i.e., u = — tt^^tttS 

° ' ' voles'") 

where vol(§ ) is the surface area of S . 

Lemma 7.3. Let d € N and let B be a matrix with entries bij £ M, i,j = 1, . . . , d + 1. 
Then 

^ ^ d+l 

j {Bv,v)^^^d.{v) = j^^Y.bn- 

Hence for zi,Z2 G M^+i we get J^a izi,v)^,^^{z2,v)^^^ diy{v) = ^ {zi,Z2)^^^. 

Proof. Define X : M"'+i -^ M'^+i as X{^) = B(,(e M"'+^ The Gaussian integral formula 
then implies 

/ {X{v),v\^^dS{v)= f div(X)(0de = voWi(i?i)^6,,-. 

Here vold+i{Bi) denotes the Lebesgue volume of the {d + l)-dimensional unit ball Bi. By 
using the well-known relation vo^S"') = {d-\- 1) volrf+i(i?i) the first claim follows. For the 
second statement just set B = Z2 zf . D 

Finally, in order to do numerical simulations we have to compute some of our underlying 
objects in local coordinates. Therfore, recall the following spherical coordinate system 
{U,x) given by x = tJ^ and U = Im(rd) with ri(6'i) := (cos(6ii) sin(6'i))'^, 0i G (0,27r), 
and inductively for d G N, d>2, 

^"-y cos(^,) j,^dG(0,vr). 

Thus the Riemannian metric on U is determined by g^- = {dgTci,dg.Tii) ■ So the density 

-' V * ^ / cue 

of the Riemannian volume measure y^ with g = det ( (^l,- ) ) , denoted by Qd, is given in 

this coordinate system by Qd = Ylj \do^Td\ since ^j- = for i / j. This yields the formula 
Qd = Qd-isin'^~^{9d), d>2. In particular, 
d 

-Qf.Bd = U - 1) cot{ej)Qd, den, j = i,...,d. 
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Define n 



'3 •" 



n 



id) ._ 



7. 1 , The gradient 



l^e^Trfl" dg^Td, j = 1, . . . ,d. Let V be as in Lemma 
of V is computed in local coordinates as gradgdF = Z]i,7=i 5*''^ ^ where (g'*-') is the 
inverse matrix of (gij). Thus 






9^ 



(7.27) 



Here Q 



3 ■~ 



\de,Td\ 



n«=j+i sinie ) ■^lisre the empty product in case j = d is defined to 



be equal to 1. Finally, for Agd = ^ Ef,j=i ^ 



^g^Jgf-) it holds 



A 



Sd 



j=l 3 



E^I 









(7.28) 
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